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Abstract 

We solve, using localization, for the large- N master field of M = 2* super- 
Yang-Mills theory. From that we calculate expectation values of large 
Wilson loops and the free energy on the four-sphere. At weak coupling, 
these observables only receive non-perturbative contributions. The an- 
alytic solution holds for a finite range of the 't Hooft coupling and ter- 
minates at the point of a large-iV phase transition. We find evidence 
that as the coupling is further increased the theory undergoes an infinite 
sequence of similar transitions that accumulate at infinity. 



*Also at ITEP, Moscow, Russia 



1 Introduction 



Many important insights have been obtained from considering gauge theory dynamics in the 
large- N limit. Perhaps the most spectacular one is the direct relationship to string theory 
through the holographic duality. In addition to the large- N limit, holography also assumes 
that the coupling constant is large. This is not necessary for the duality to hold, but the 
holographic description is simple (and therefore useful) only in the strong-coupling regime. 
Needless to say, any direct computations in a strongly coupled field theory are difficult, but in 
supersymmetric Yang-Mills theories many non-perturbative results became available with the 
development of powerful non-perturbative methods, such as integrability and localization. 

In this paper we will study J\f = 2* theory, where the localization methods can be readily 
applied [TJ. The M = 2* theory is the unique massive deformation of M = 4 super- Yang- 
Mills (SYM) that preserved J\f -2 supersymmetry. Localization reduces the partition function 
of N = 2* SYM to a finite-dimensional matrix integral [lj. Solving the M = 2* localization 
matrix model at large N and large 't Hooft coupling opens an avenue for confronting direct 
field-theory calculations with holography in this non-conformal setting, since the supergravity 
background dual to M = 2* SYM is also explicitly known [21 E]. The strong-coupling solution 
of the matrix model, obtained in [4), appears to be in the full agreement with the predictions 
of the holographic duality. 

Our goal here is to investigate the M = 2* theory away from the supergravity limit, at 
arbitrary values of the 't Hooft coupling. We will solve the localization matrix model non- 
perturbatively by a method proposed by Hoppe [U E]- Quite surprisingly, the solution thus 
obtained behaves well if the coupling is not very big. At some finite value of the coupling the 
solution terminates, signaling a transition to a new phase. Phase transitions of this type are 
ubiquitous in the large- N theories [7], and in the early studies of the AdS/CFT correspondence, 
it was even hypothesized that the weak-coupling (perturbative) phase of M = 4 SYM would 
be separated from the strong-coupling (string) phase by a large- N phase transition [Sj. This 
has never been confirmed, and it is now quite firmly established that the coupling constant 
dependence in M = 4 SYM is smooth all the way from weak to strong coupling. We will see 
that in the M = 2* theory this is no longer the case. There is a phase transition separating strong 
and weak coupling regimes. Preliminary numerical investigation of the strong-coupling phase 
indicates a very complex structure, with infinitely many phase transitions accumulating at 
infinite coupling. In the asymptotic strong-coupling limit one nevertheless arrives at the simple 
Wigner-type solution found in [4], that agrees precisely with the supergravity predictions [3j. 

2 Af=2* theory 

The M = 2* theory is a relevant perturbation of maximally supersymmetric M = 4 SYM by a 
combination of dimension two and dimension three operators that preserves half of the super- 
symmetry. The field content of M = 2* SYM therefore is the same as in the M = 4 case: there 
are six scalars $i, . . . , $4, $, $' and four Majorana fermions in addition to the gauge field A^. 
All fields are in the adjoint of the gauge group, which we take to be SU(N). The relevant 
perturbation adds equal masses to $/ and their superpartners, which together form an Af = 2 
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hypermultiplet. The vector multiplet contains $, $' and the gauge field, and remains massless. 
We denote the mass scale of the Af = 2* theory by M. 

In the IR, at energies much below M, the hypermultiplet decouples and the theory flows to 
the pure Af = 2 gauge theory. One can thus view Af - 2* theory as a particular UV regularization 
of Af = 2 SYM. The hypermultiplet mass M then plays the role of a UV cutoff. The pure Af = 2 
theory, unlike Af = 2*, has a non-zero beta functions with an associated dynamically generated 
scale 

A=e-^M, (2.1) 

where A = gy M N is the Af = 2* 't Hooft coupling. It is important to emphasize here that the 
flow picture makes sense only if the UV and IR scales are widely separated (A « M), which 
requires the 't Hooft coupling to be small. 

One can improve the low-energy effective field theory, that arises upon integrating out 
hypermultiplets, by including operators of higher dimensions. An operator of dimension d 
will enter the effective action with a potentially calculable coefficient proportional to M 4 ~ d . 
Its contribution to observables will thus be suppressed by (A/M) d ~ 4 . These standard OPE 
arguments suggest that at weak coupling any observable A in Af = 2* SYM has a double 
expansion 

o oo 9 

A = M A e~— Y,C n e~—, (2.2) 

n=0 

where A is the dimension of A, and the OPE coefficients C n themselves are power series in 
A. We have taken into account that operators that enter the effective action have even integer 
dimensions. 

Using localization we will compute a number of observables non-perturbatively. We will see 



that their weak-coupling expansion has precisely the expected OPE structure (2.2), with OPE 
coefficients calculable to all orders in A 2 /M 2 . In the cases we consider the OPE coefficients 
C n turn out to be numbers that do not dependent on A, perhaps due to supersymmetric 
non-renormalization theorems, as only supersymmetric observables can be computed using 
localization. 

The mass deformation of Af = 4 SYM triggers symmetry breaking of SU(N) to [/(l)^ -1 , 
pushing the theory into the Coulomb phase. The symmetry breaking is characterized by a 
diagonal vev of the scalar in the vector multipletj^] 

($) = diag (ai, . . . , a N ) . (2.3) 

The distribution of symmetry-breaking eigenvalues plays the role of the large- N master field 
[10J of M = 2* theory. In the large-A^ limit the distribution of eigenvalues is described by a 
continuous density 

1 N 

P^) = irZHx-a l ), (2.4) 

ly 8=1 



1 More precisely, different vacua are characterized by a diagonal vev of the complex field $ + The 
low-energy dynamics on the moduli space of vacua is described by the Seiberg-Witten theory [5]. Since the 
localization matrix integral of pQ is convergent when the contour of integration runs along the real axis, we shall 
assume that all a,'s are real. 
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which does not fluctuate and is determined by the classical equations of motion. 

As usual in the large- N matrix models, the eigenvalue density p(x) is different from zero 
only on a finite interval [-p, p]. The largest possible eigenvalue p is an important dynamical 
quantity that sets the scale of the symmetry breaking. For instance, the heaviest W-boson 
mass is equal to 2p. As we shall argue shortly, p also determines the behavior of large Wilson 
loops. 

The Wilson loop expectation value is defined as 



W{C) 




J ds (ix^A^ + 

c 



(2.5) 



where the coupling to the scalar is dictated by supersymmetry. In a first crude approximation 
(which nevertheless turns out to be exact for sufficiently large contours), we can just substitute 



the master field (2.3) for $ and neglect all quantum corrections: 

l N r 

W(C) = —Te aiL = / dxp(x)e xL ^ e^ L . (2.6) 
N i=1 i 

Here L is the length of the contour C and the last equality holds as far as L » fjr 1 . If the 
theory is compactified on S 4 , and the loop runs along the big circle of the sphere, the first step 
in this calculation, replacement of $ by its vev, is not an approximation. The circular Wilson 
loop is a supersymmetric observable and can be computed by localization. Its vev is thus given 
by the exponential average in the resulting matrix model pp. In the large- N limit, the quantum 
average over the ensemble of dj's is replaced by the classical average over the density p(x). If 
we further take the decompactification limit L -> oo, the integral over the eigenvalue density is 



peaked at the upper bound of integration and the estimate in (2.6) gives the Wilson loop vev 
with an exponential accuracy. The logarithm of the circular Wilson loop thus grows linearly 
with the size of the circle. Since the behavior of large Wilson loops should be fairly universal 
and largely independent of the contour's shape^J we conclude that Wilson loops in Af - 2* SYM 
satisfy perimeter law with the coefficient of proportionality between \nW(C) and L equal to 
p. 

On dimensional grounds, the symmetry breaking scale p should be proportional to M, the 
only mass scale in the problem: 

p = g(X)M. (2.7) 

Our main goal will be to determine the eigenvalue density p(x) and the coefficient of propor- 
tionality g(X). 

The eigenvalue density p(x) and the symmetry breaking scale p are known in the two 
limiting cases of very large and very small A. At weak coupling, the only relevant scale is A 



from (2.1), and the dependence of g on the 't Hooft coupling is dictated by the asymptotic 



freedom of the low-energy Af = 2 SYM: 

#(A)=2e^ (A-»0). (2.* 



2 At strong coupling this can be explicitly verified using the holographic description 
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The coefficient of proportionality was computed in [TT]. The eigenvalue density obeys the 
inverse square root law at weak coupling: 

P(s) = X ~, ( 2 - 9 ) 

%\J [l l - X 1 

This result was derived in [12] from the Seiberg-Witten theory, and can be also obtained from 
localization 

The eigenvalue density at strong coupling is the Wigner semicircular distribution: 

p{x) = —J^-xi (A-oo), (2.10) 

TTfI z 

with 

g(X) = ^ (A-oo). (2.11) 

The strong-coupling results follow from the probe analysis of the supergravity dual [3J, or can 
be derived directly from the localization matrix model in full agreement with the predictions 
of holography [I] . 

We will analyze the localization matrix model in the intermediate regime of finite 't Hooft 
coupling. The first step in the localization approach of [TJ consists in putting the theory on the 
four-sphere of radius R. Here we regard this step merely as an IR regularization and will be 
mainly interested in the decompactification limit R -> oo. Nevertheless, considering the theory 
on a sphere introduces an extra observable, the free energy, which we can also calculate from 
the matrix model. Since the vacuum energy is zero by supersymmetry, the free energy is not 
an extensive quantity. It turns out that at large radius the free energy has the form 

F = -InZ = -N 2 M 2 R 2 (In MR + f(X)) + (R°) . (2.12) 

The leading R? In R piece is of purely one-loop origin. The coefficient in front is not renormalized 
by quantum correction^} The A-dependence starts at 0(R?) and is parameterized by a single 
function of the 't Hooft coupling /(A). This function is known in the two limiting cases [TT] 14"]: 



5 In A + 2f + 7 - m47r + . . . (A — oo) , 



where 7 is the Euler constant. 



3 Saddle-point equations 

The master field of M = 2* theory is determined by the saddle-point equations of the local- 
ization matrix model. These equations greatly simplify in the decompactification limit, and 
can actually be derived from simple heuristic arguments. We will later re-derive them more 
rigorously directly from the localization integral. 

3 It would be very interesting to check this non-renormalization property by computing the free energy at 
strong coupling from the on-shell action of the supergravity solution, which unfortunately is not known for the 
theory compactified on S 4 . 
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3.1 Heuristic derivation 



The mass spectrum of M = 2* SYM in the background (2.3) is given by 

mjj = \a,i-aj\ (3.1) 
mjj = \a i -a j ±M\, (3.2) 

where m v and m h are masses of vector and hypermultiplets. At weak coupling, when fi « M, 
the vector multiplet masses are of order m v ~ fi, while hypermultiplets are heavier: m h « M. 
However, the ratio g = [i/M grows with A, and at strong coupling the ordering of scales is 
reversed: // » M. Then all the masses are determined by the largest scale: m v ,m h ~ //. 

These estimates apply to typical states. In addition, the spectrum contains a number of 
exceptionally light states with parametrically small masses. At weak coupling these are the 
W-bosons from the vector multiplet with \i - j\ ~ 1 (as opposed to \i - j\ ~ N). The masses of 
those states scale as 1/N in the large- A^ limitj^J mjj ~ 1/N. 

Interestingly, hypermultiplets can also be parametrically light as soon as fi becomes suffi- 
ciently large. If the width of the eigenvalue distribution exceeds M (2/i > M), there will be 
eigenvalues a; and aj with a^-aj « ±M, for which m^. ~1/N. Resonance phenomena associated 
with these light states have dramatic consequences for the structure of the master field. As we 
shall see, the appearance of the first resonance leads to a large- A^ phase transition. The critical 
coupling of the transition is determined by the condition 2/i = M, i.e. 

g{K) = \- (3.3) 

Appearance of new resonances leads to secondary phase transitions. This happens each time 
2(1 crosses an integer multiple of M: 

t{W)-\- (3-4) 

We will explicitly see the first transition in the analytic solution of the matrix model, and we 
will provide numerical evidence for the secondary transitions. 

To derive the saddle-point equations for the master field, we need to minimize the effective 



action for the background (2.3) on the sphere of large radius. The rigorous justification of 
this procedure comes from the analysis of the localization matrix model of [T]. But since the 
effective action is one-loop exact, it is easy to give a heuristic derivation based on the form of 



the spectrum (3.1), (3.2) and simple physical arguments. The effective action consists of the 
classical piece, originating from the supersymmetic coupling of the scalar $ to the curvature 
ofS 4 : 

S cl =-\- f d^x^glZtr® 2 , (3.5) 



4 As explained in [T2], the masses of the lightest states actually scale as 1/N 2 . This happens because the 
density blows up at the edges of the eigenvalue distribution. However, the entropy associated with the light 
states is small in the large-TV limit - there are O(l) states with m ~ 1/N 2 , O(N) states with m ~ 1/N and 
0(N 2 ) states with m ~ 1. 

5 The supersymmetric coupling differs by a factor of 3/2 from the conformal 7Z/6 coupling, as follows e.g. 
from eqs. (2.1), (2.3) in [T3]. 
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and the one-loop corrections from integrating out massive fields. The contribution of a field of 
mass m must be proportional to m 2 R? In m 2 R 2 . The coefficients of the vector and hypermulti- 
plet contributions can be fixed by requiring that (i) at M » |aj| the only effect of the one- loop 
corrections is the coupling constant renormalization with the correct beta function of the M = 2 
theory, and (ii) at M = the one-loop corrections vanish. 

Taking into account that the volume and the scalar curvature of the four-sphere are equal 
to Vo^S* 4 ) = 87r 2 i? 4 /3 and 1Z = 12/R 2 , the classical part of the action becomes 

8tt 2 /V n 

Sci = ^£a^ 2 - (3.6) 

A i=i 

The one- loop quantum corrections are uniquely fixed by the conditions (i) and (ii), resulting in 
the following effective action: 



S cS [a] 87r 2 N ^ 2 ^ 



- (en - aj + M) 2 In (ai - aj + M) 2 R 2 



+ - (ai - aj - M) 2 In { ai - aj - M) 2 R 2 -- (a* - aj) 2 In (a, - ajf R 2 ] . (3.7) 
4 2 J 

The i?-dependence on the right-hand side factors out, in the sense that the independent term 
does not depend on a^. This leads to the one-loop exact R 2 lni? term in the free energy discussed 
earlier. 

The full free energy is then obtained by minimizing the effective action in af 

1 dS cS 8vr 2 1 

= ^7 = X fll - n g G{a * - a ^ (3 ' 8) 

where 

G(x) = (x + M)ln\x + M\ + (x-M)ki\x-M\ -xlnx 2 . (3.9) 
This set of equations is equivalent to an integral equation on the eigenvalue density: 



^ o 9 



dyp(y)G(x-y) = —x. (3.10) 



The integral equation greatly simplifies when A is large, because then M « fi, and the kernel 
can be Taylor expanded in powers of M: G(x) = M 2 /x + . . . In this approximation, the equation 
is equivalent to the saddle-point equation for the Gaussian matrix model with an effective 



coupling A e ff = AM 2 , and is solved by the Wigner distribution (2.10) with \i given by (2.11) [1]. 
At weak coupling, the opposite inequality holds: M » /i. The first two terms in the kernel 
are then approximated by 2x In M and can be transferred to the right-hand-side, leading to the 
coupling constant renormalization: 8tt 2 /X -> 8tt 2 /X - 21nM. The resulting integral equation is 



solved by (2.9), (2.8) [11]. We will solve the saddle-point equation exactly in a finite range of 



A in the next section. 
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We will also compute the function /(A) defined in (2.12). In order to do that, we can 
differentiate the effective action in A: 

/,(A) = ~WWr 2 d W = ^W^^ = ^pf dxp{x)x2 ~ iw 2 ^ ■ (3 - n) 



The mass dependence can be eliminated from the above equations by rescaling x -»■ xM, 
H -> gM. It is thus possible to set M to one by choosing appropriate units. However, it will 
be convenient to keep it as a parameter, because the equations are easier to solve for complex 
values of M. The solution for real M can be then obtained by analytic continuation. 

3.2 Localization 

A more rigorous derivation of the saddle-point equations starts with the localization partition 
function of M = 2* theory pQ: 



d a 11 rr ( L rr) 6 1 ' 3 ' 12 

H(di - cij - M)H{ai - dj + M) 



where 

%)=niv e ^ ( 3 - 13 ) 

n=i v n / 

The instanton partition function Zi nst is a known [14], albeit complicated function of the eigen- 
values Oj, the gauge coupling and the number of colors. Fortunately, instantons are suppressed 
at large iV and the instanton factor |J?i nst | 2 can be set to one as soon as we are interested in 
the large-iV limit. In principle this step requires justification, which will be given elsewhere 
[TT| [15] . Here we just assume that the instanton factor can be dropped. This leads to a major 
simplification of the localization integral, which at large N can be computed by the saddle-point 
method. 

The saddle point equations for the eigenvalue density read 
n 

j d VP(y) (^^-K(x-y) + ^K(x-y + M) + ±K(x-y-M)} = 8 ^x, (3.14) 

-A* 

where 



, R 2 x - y 

H'(x) 



These equations are valid for any R and any A and were analyzed in different corners of the 
parameter space in [HI [HI H] . 

We are going to concentrate on the decompactification limit R -*■ oo. In this limit the 



Hilbert part of the kernel is suppressed by the 1/R 2 factor, and the summation over n in (3.15) 
can be replaced by integration over v = n/R: 

oo 

K(x)^2x f dv{ ^—Ac±2x\rR\x\. (3.16) 

J \ V V 1 + X 2 ) 

1/R 
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Here R can be regarded as an IR cutoff that regularizes the contribution of light W-bosons 
to the effective action. The cutoff dependence however disappears in the integral equation, 
reproducing (3.10), (3.9). 



4 Solving the matrix model 



In order to solve the saddle-point equation (3.10), we differentiate it twice. After the first 
differentiation we get: 

(x - y) 

Differentiating again, we arrive at a singular integral equation: 



/ dyp(y)\n y J> =—. (4.1) 



dyp{y) -0, (4.2) 



(x - y)[(x - yf - M 2 ] 



that can also be written as 



-fdyp(y)l— 1 —— ^rrUo. (4.3) 

J yHKyj \x-y x-y + M x-y-M) y ' 

-v- 



The general theory of integral equations with a singular kernel [T7] tells us that (4.2) has a 



one-parametric family of solutions, which behave at the boundaries of the integration domain 

as 

, . const , , t s 

p{x) ~ — = (x -+ ±n) (4.4) 

y/H T X 

In this case the free parameter is easy to identify, since the equation is invariant under p(x) -*■ 
const p(x). The properly normalized solution is consequently unique and exists for any p. The 



integrated form of the saddle-point equation, eq. (4.1), then determines p as a function of A. 



We will solve the integral equations (4.2 ), (4.1 ) by analytic continuation in M. The integrand 



in (4.1 ), as a function of M, has a logarithmic cut between ±(x-y), and the analytic continuation 
is unambiguous only if M 2 > (x - y) 2 throughout the whole region of integration. This requires 
M to be bigger than 2p. This is the same condition that we encountered in the discussion 



of resonances in the beginning of sec. 3.1 In fact, the branch points at M = ±(x - y) are 
directly related to the resonance due to a nearly massless hypermultiplet. The appearance of 
this resonance at p = M/2 leads to a failure in the analytic continuation and, as we shall see, 
makes the dependence of p on A non-analytic, causing a transition to a new phase. 



4.1 Solution of the auxiliary problem 

For M 2 > (x - y) 2 the transformation M ->■ iM changes the imaginary part of the logarithm in 
(4.1 ) by 7r, which has to be compensated by an imaginary shift of the coupling on the right-hand 
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side. The analytic continuation thus boils down to the following substitution: 



M iM, —— = —^ + i-k. 
A A 



(4.5) 



After the analytic continuation, the saddle-point equations (4.1), (4.2) take the form 



J dyp(y) 

I- 1 



In 



M 2 + (x-yy 8tt 2 



(x-yY 

dyp(y) 



A 



/ (x-y)[Mi + (x-y) 2 ] 



0. 



(4.6) 
(4.7) 



We first solve them assuming that M and A are real, and then analytically continue back to 
real M and A. 



An integral equation with the same kernel as in (4.7), but with the right-hand side linear in 
x was solved in [5j E] . The method of (5J E] can be readily adapted to the case at hanc^ 
idea consists in considering a generalized resolvent 



The 



G(z) = f 



dyp{y) 



(4- 



This function is analytic in the complex z plane with two cuts from -p±iM/2 to p±iM/2. In 
the vicinity of the upper cut: 




A* 



dyp{y) 



{x-y)+M 2 



dyp(y) 



M 



J (x 



y)[M* + (x-y) 2 ] 



(4.9) 



The last term vanishes in virtue of the integral equation (4.7), which means that G(z) is real on 
the two sides of the cut. Since the reality condition is fully equivalent to the integral equation, 
it suffices to construct a function which has the same analyticity properties as G(z) and which 
is real on the same contour. The discontinuity of this function across the cut automatically 
satisfies the integral equation. 

The locus on which G(z) is real forms a closed curve that connects +ioo with +oo along the 
contour shown in fig. [TJ Instead of G{z) it is actually more convenient to deal with 1/G(z), 
which monotonically increases from -oo to +oo along this contour. Moreover, the imaginary 
part of 1/G(z) is positive in the upper-right quadrant of the complex plane with the cut and 
consequently 1/G(z) maps the interior of the contour shown in fig. [ljonto the upper half-plane. 
Such map is unique, and its inverse can be explicitly constructed in terms of elliptic integrals. 
This is the idea of the method of [3 [6]. 



6 It is impossible to just set the force term to zero. The solution obtained in [SJ [3] is singular in this limit. 
To get the right result one should properly take into account the difference in the boundary conditions. 
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-00 




+00 



Figure 1: The function G{z) is real on the contour shown in the solid line. The arrows correspond to 
the positive gradient of 1/G(z). Also displayed are the values of 1/G(z) at the turning points on the 



curve. The dashed line is the contour of integration in eq. (4.27). 



To construct the map from the upper half-plane to the domain shown in fig. [TJ we need 
to know the values of the function l/G at the turning points of the curve. At infinity, G(z) 
behaves as 1/z 2 , as follows from its definition. For later purposes we will need few more terms 
in the expansion: 



G{z) = 7* + i [ 3{x2) - x) + h (if - + 5(x4) 1 + •• ( ^ oo) • (410) 



Since the density blows up at the endpoints of the eigenvalue distribution, according to (4.4), 
eq. (4.9) implies that 

1 1 / 

-0, —— ^ = +0. (4.11) 



G^f + n + ie) ' G{if+n-ie\ 
Finally, we denote the values of l/G at three other turning points by -£ 3 , £ 2 and £i: 

It can be shown that > and that £i > £2- Otherwise, £j's are unknowns, to be determined 
later. 
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With these data at hand, the map z -> G 1 (z) is given by an elliptic integral: 

i 

dpp 



- f 

2 I V(p + 6)(p-6)(p-6) 



(4.13) 



The constant in fro nt (1/ 2) ensures that 2dz = \fGdG~ 1 at infinity (z -> oo), in accord with the 
large- z asymptotic (4.10). The correct behavior near the endpoints, G(z) ~ (iM/2 + fi- z)' 1 / 2 , 
emerges here by construction, since dz ~ dG~ l G~ l when G' 1 -»■ 0. 

To find the parameters of the map, we integrate dz along all the finite segments of the 
curve on which G(z) is real and compare the results with the geometric data from fig. [Tj 



SI 



dpp 



6 



M 



«2 



dpp 



/v(p+6)(6-p)(6-p) 



2/x. 



(4.14) 
(4.15) 
(4.16) 



These equations fix 6 in terms of M and \i and together with (4.13) completely determine the 
solution of the integral equation (4.7). As expected, the solution exists for any fi. Since the 
normalization condition on the density is built in the construction, the solution is also unique. 
Let us rewrite the conditions that determine 6 in terms of the canonical elliptic integrals: 



2 (6 + 6) #'-26^' = MVfi + fa 

(6 + 6)#-6# = 

(ei+e 3 )^)-ei^) = Mi +6 ■ 



(4.17) 
(4.18) 
(4.19) 



Here E = E(m), K = K(m), E(ip) = E(ip\m), F(<p) = F(if\m), E' = S(l-m), and K' = K(l-m) 
are elliptic integrals of the first and second kind with the modulus and the angle given by 



m 



sin 2 ip 



6 + 6 
6+6 
6 

6 + 6 



(4.20) 
(4.21) 



We can express 6 and \i as functions of the modulus parameter by inverting these equations. 
Using the Legendre's identity, 



KE' + EK' - KK' 



71 



(4.22) 
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we get for the £j's 





M 2 




6 = 




KE 


7T 2 






M 2 


K [E 


6 = 






T 2 " 




M 2 


K (K 


6 = 


T 2 " 



(l-m)K] 
-E), 



From (4.19) we then have: 



M 



p = -(KE( V )-EF( V )) 



71 



where the angle is determined by (4.21): 



sin 2 cp 



K-E 
mK 



(4.23) 



(4.24) 



(4.25) 



Varying m we get the solutions of the integral equation for all possible p. 

We want, however, to find p as a function of A. In order to do that we should impose the 
integral condition (4.6). This can be done as follows. Let us integrate the resolvent from some 



point on the upper cut to infinity. From the definition (|4.8|) we get 

■kJ"""^ 



I 



dzG(z) 



i In 



x-y 



y + iM 

%M 2iM 



L- fdy P (y)ln ^ + (x -f . (4.26) 



(x-yY 

For real x the first term on the right-hand side is real and the second term is imaginary. 
Comparing to (4.6), we see that the integral form of the saddle-point equation is equivalent to 
the condition 

r An 2 

/ dzG(z)=-^r, (4.27) 
J. AM 



Im 



which should hold for any real x between and \i. The contour of integration here is arbitrary. 
For instance, we can integrate along the dashed line in fig. [T] On the horizontal segments of 
this line both dz and G(z) are real, and these segments do not contribute to the imaginary 
part of the integral. This, in particular, insures that the result does not depend on x. 

The integral in the last equation can be computed with the help of the differential form of 
the map (4.13): 



dz 



dpp 



2v / (p + £ 3 )(p-6)(p-6) 

Multiplying by G and integrating along the interval of p ranging from £ 2 to £i we get: 



(4.28) 



?2 



dp 



8tt 2 



(4.29) 
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Written in terms of the standard elliptic integrals this equation becomes, upon using (4.23), 



K' 

~K 



An 

X 



(4.30) 



The three equations above, (4.24), (4.25) and (4.30), give \x as a function of A in a parametric 



form. It is possible to solve for li directly by expressing elliptic integrals in terms of theta 
functions. The definitions and some properties of the theta-functions are listed in appendix |A} 
The theta-functions are defined on a torus with the modular parameter r = iK'/K. This is very 



convenient, since in virtue of (4.30) the modular parameter is directly related to the coupling 
constant: 

(4.31) 



Q 



_ 47T 

e x 



Using the formulas of appendix [A] the parameters of the solution can be written as explicit 
functions of the coupling: 



6 
6 



M 2 

rT 
~12~ 

M 2 
12 



(E 2 -6 4 2 + 26l) 
(E 2 + 26l-6l) 

(-E 2 + et + el). 



We use the following shorthand notations: 9 a (z) = 9 a (z\q), 6 a 
Finally, the elliptic modulus and the angle are given by 



a ((%) and E a 



ni 



sin 2 Lf 



9J 

291 -Q\-E 2 
391 



(4.32) 
E a (q 2 )- 

(4.33) 
(4.34) 



Another form o 



The equation (4.24) then gives ii directly as a function of A. 

the equations for fi is obtained by replacing the elliptic integrals in 



(4.24) 



(4.25) by the theta-functions: 



goo 

0\{v) 



M d 
~2 ~di 



ln9 4 (v) 
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01 



E 2 



39 2 9 2 



(4.35) 
(4.36) 



These two equations again determine fi in the parametric form. The last equation can be 

TrF(ip) 



inverted by using (A.13): 



2K 



F(<P) 



(4.37) 
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In addition to the largest eigenvalue //, it is also possible to compute the averages (x 2n ). We 



proceed as follows: equation (|4.13|), written in the differential form, becomes 

I 



dz = - 



dG 



(4.38) 



2G3/2 ^(1 + 60(1-60(1-60) ' 

The expansion of this last expression in powers of G can be compared, order by order, with 
the expansion of the resolvent at large z, given by (4.10). This gives a set of algebraic relations 
between (x 2n ) and that can be used to recursively solve for (x 2n ). Expanding to the first 
two orders, for example, we get: 



(x 2 ) = _L(M 2 -46-46 + 46) 



(4.39) 



(x4) = 240 ^ ' 40A ^ (6 + 6 ~ 6) + 56 (6 + 6 " Cs) 2 " 8 ^ + ^ + ^ • (440) 



Substituting the explicit form of 6's from (4.32), and using the theta-function identities listed 
in appendix [X we get (x 2 ) and (a; 4 ) as explicit functions of A: 



, *i M 2 , v 
(^ 2 ) = ^-(l-^) 

M 4 

<x 4 ) = — (l0£ 2 2 - 30£ 2 -£ 4 + 2l). 



(4.41) 
(4.42) 



We will use the first of these equations to calculate the free energy in the next subsection. It 
is straightforward to push the expansion in principle to any desirable order and express (x 2n ) 
as a combination of Eisenstein series of up to the 2n-th degree. 



4.2 Free energy 

In solving the auxiliary problem we assumed that M and A (and consequently q) are real. The 



solution of eq. (4.36) for v then is also real giving a real-valued /i, as expected. Upon the 



analytic continuation to real A, A becomes complex and q becomes pure imaginary: 



q = ie 



X 



(4.43) 



It is not immediately clear from the equations defining /i that it will remain real-valued after 
the analytic continuation. We will return to the reality condition on /x shortly. For now on 
we are going to discuss the analytic continuation of the free energy, which is more or less 
straightforward. 



The derivative of the free energy in the coupling is given by (3.11). Taking (x 2 ) from (4.41 ) 



we get, upon the analytic continuation (4.5) 



f(A) 



2tt 2 
3A 



2 (E 2 (-p)-l) 



(4.44) 
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where we have introduced 

p = -q 2 = e~^~ . (4.45) 



This equation can be integrated with the help of (A.5), yielding 



3A 12 

where r](x) is the Dedekind function. Explicitly. 



/(A) = -^-^ + 21nij(-p) , (4.46) 



/(A) = 2J> l-C-lfe"^ . (4.47) 

n=l V / 

This surprisingly simple formula is one of our main results. Few remarks about this formula 
are in order. First, it nicely matches with general expectations for the weak-coupling expansion 



(2.2). To the first few orders, 



f , n „ 16n 2 8 24tt 2 7 327T 2 12 _ 40n^ 48tt 2 

/(A) = 2e A-3e a +- e a — e * + — e a -4 e a +... (4.48) 

3 2 5 



The first term agrees with the weak-coupling calculations of [llj . quoted in eq. (2.13). The 
rest can be identified with higher orders of the OPE. Interestingly, only non-perturbative terms 
appear in the expansion, and all the OPE coefficients turn out to be simple rational numbers. 
This is perhaps a consequence of supersymmetry non-renormalization theorems. 

It is tempting to continue the free energy to the strong coupling and compare the strong- 
coupling expansion with expectations from string theory. Our subsequent results, however, 
show that this calculation is meaningless because of the phase transition that happens before 
one reaches the strong-coupling regime. The free energy itself does not contain any signs of 
non-analyticity in A, which is typically the case for discontinuous phase transitions, but it is 
possible to see that something wrong happens at sufficiently big A even at this stage. The 
derivative of /(A) is proportional to the average (x 2 ) and consequently must be positive. But 



the right-hand-side of eq. (4.44) becomes negative for A > 88.8, signaling an instability of the 



weak-coupling phasd^j The reason for this unphysical behavior is violation of positivity of the 



density p(x). At the critical point determined by (3.3), the density becomes complex, such 
that its moments are no longer positive-definite. 

Let us finish with few speculative comments on the exact formula for the free energy. The 
function /(A) coincides with the free energy of a free 2d boson with twisted boundary conditions, 
and the temperature related to the 't Hooft coupling. We do not know if this observation has 
any physical significance. The result is also strikingly similar to the free energy of a toy model 
considered by Pestun [1] , in which the e-parameters used in the localization of the path integral 
for M = 2* are adjusted to cancel the perturbative contribution. The resummed instanton 
contributions then lead to a similar formula for the free energy, with gyyi playing the role of A 
and a different proportionality constant (N instead of N 2 M 2 R 2 ). 



7 For higher moments of the density, the positivity is violated even earlier: (a; 4 ) becomes negative at A « 61.2, 
(a: 6 } becomes negative at A « 53.8. As we shall see, the phase transition happens as A c « 35.4. 
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q c n 



1 

Figure 2: Analytic continuation contour in the g-plane: the solution of the auxiliary problem with real 
A corresponds to the dashed line. The physical domain of real A is shown in solid line. The dotted line 
corresponds to the unphysical solution obtained by analytic continuation past the phase transition. 



4.3 Analytic continuation 



(4.36) 



We n o w proc e ed wi th the analytic continuation of eqs. (4.35) 
(4.37), (4.33), (4.34), to real A and M. To do so we start with large real A (q 



or equivalently (|4.35|) 
1 



decrease 

q to zero and then make it pure imaginary. The continuation path in the complex q plane is 
shown in fig. [2] At q = 0, we find from eqs. (4.33), (4.34) and (4.37) that m = and v = cp = 7r/4. 



Eq. (4.35) then gives 



/i~2Mg = 2Me~^, (4.49) 

which is real for physical values of parameters and reproduces the weak-coupling result ( |2.8 ) 
[TT] for the function g(X) = /i(A)/M. It is straightforward to compute few more orders by 
expanding in higher powers of q: 



g(\) = 2e 



4irf 
A 



4e ~ 



16e 



58e * 



44n z 

324e- — 



1856e~ 



+ ... 



(4.50) 



The weak-coupling expansion again has the form (2.2) expected from the OPE. The OPE 



coefficients appear to be numbers, independent of the coupling, probably also due to supersym- 
metry non-renormalization theorems. Following the discussion of section 2, the circular Wilson 
loop will be given b}|^] 

W(C) * e 2nMRg W . (4.51) 

The g-series for \x contains only odd powers of q and, since \i is linear in M, it is indeed 
real after analytic continuation to imaginary q and imaginary M. This can be demonstrated 
directly from (4.35), (4.36). First we notice that to the linear order in q, v ^ ir/A-q. Observing 



that q is pure imaginary after the analytic continuation, we can look for solutions of (4.36) of 
the form 

v = - + - } (4.52) 
4 2' v ' 



3 We expect this behavior to hold for any large contour. 
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Figure 3: The path of the analytic continuation in the u-plane. In order to make the horizontal lines 
straight in the last plot the imaginary axis must be rescaled by 1/A. 



where s is real and negative (fig. [3]). 

To get rid of the imaginary i in q, it is convenient to perform the Landen transformation 
v -> 2v, q -*■ q 2 in the theta-f unctions, followed by a modular transformation q 2 -*■ -q 2 , and a 
shift by a half-period 2v -> 2v + tt/2. The transformation rules for the theta-functions can be 
found in |18j . After this chain of transformations (4.35), (4.36) become: 

1 d 

9 



2 ds 



ln(0 4 (O)0 3 («s) + «02(O)0i(is)) 



9 2 (Q)9 4 (0)9 2 (is) - ielMOxiisfaiis) _ 2 4 (O) - 6f(0) + E 2 



eiiis) 



602(0)04(0) 



(4.53) 
(4.54) 



where now 6 a (z) = 9 a (z\p), E 2 = E 2 (-p) with p denned in (4.45). We also wrote the equations 
in terms of the dimensionless ratio g = \i\M. These equations look more complicated than 



(4.35), (4.36), but now the modulus parameter of the theta-functions is real. Since 



i9i(is),9 2 (is),9 3 (is),94(is) e R for s € R 



(4.55) 



eq. (4.54) is a real equation for a real variable s. The function g(X) then is also manifestly real. 

The contour of analytic continuation in the f-plane, consequently, has the shape shown in 
fig. [3j v is real and changes from to 7r/4 for the auxiliary problem. The physical regime 



corresponds to v of the form (4.52) with real negative s. The magnitude of s grows with A, 



but it is clear from the outset that it cannot grow indefinitely because the theta-functions are 
essentially periodic. What happens is that the solution for s ceases to be real for sufficiently 
large A. The critical point corresponds to the large- iV phase transition in the matrix model 
that we expect to occur on physical grounds. 
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4.3.1 Phase transition 



As discussed in section 3.1, the phase transition can be understood as a resonance phenomenon 
due to the appearance of light hypermultiplet states, with |aj - Oj| ^ M. Such states can 
only appear for A sufficiently large, because at very small A all |aj - aj\ ~ \i « M and all 
hypermultiplets are heavy. But the scale of symmetry breaking \x grows with A. The first 
resonance appears in the spectrum when \i reaches M/2. This has dramatic consequences for 
the structure of the master field leading to a large- N phase transition. 

There are many ways to see how the transition happens in the exact solution of the saddle- 
point equations described above. The easiest way to see that \i - M/2 is a critical point is to 



inspect the analytic structure of the resolvent (4.8). At weak coupling the resolvent has two 



distinct cuts in the complex plane. When fi approaches M/2 the two cuts collide, leading to 
a non-analytic dependence on A. This mechanism is quite common in large- N matrix models, 
although usually the cuts are directly related to the eigenvalue density. Here the eigenvalue 
density is concentrated on a single cut both before and after the phase transition. The cuts 
that collide are mirror images of the single physical cut and arise due to a particular structure 
of the integration kernel. 

At the technical level, the phase transition manifests itself in the divergence of G(z) at z = 
(this is the point where the cuts collide). Normally, G(0) is finite, but at the critical point it 
blows up. We actually used G(0) in solving the integral equations: the inverse of G(0) was 



denoted by £i in (4.12). The critical coupling is thus determined by the condition 



£i(A c ) = 0. 



(4.56) 



Let us demonstrate that this condition is equivalent to (3.3) and that the solution described 



above indeed hits a singularity at the critical coupling and cannot be analytically continued 
past A = A c without violating reality and positivity of the density. From the explicit expression 
for £i in terms of elliptic integrals, eq. (4.23), we see that at the critical point the elliptic E 



vanishes. Eq. (4.25) then implies that 



sin 2 <p 



at 



m 



A = A C . 



(4.57) 



But this is precisely where the elliptic integral F(cp) has a branch point. If we use (4.37) to 
solve for v, at A = A c the solution will hit a branch-point singularity, and the curve v = v(X) will 
make a turn by 90° in the complex plane (fig. |3|. Past the transition point, v will no longer 
have the form (4.52) with real s, making \x complex. 

It is also possible to see that the solution cannot be analytically continued beyond A = A c 



by inspecting eqs. (4.53), (4.54) that determine g as a function of A. The dependence of g on A 



is given there in a parametric form with the parameter v. To get the idea where the transition 



happens in the f-plane, we substitute (4.57) into (4.37) and use the standard identities for the 



elliptic integrals to compute v right at the critical point: 



v(X c ) 



vrF(arcsin^) nK (±) n i n K' tt 2n 2 i n 2nH 



2K 



2.MK 



2K 



A, 



7T 

4 



(4.58) 
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Figure 4: Image of the analytic-continuation contour in the g-plane, g = \i\M. 



where in the last two steps we used (4.30) and (4.5). 
reaches -4ii 2 /X. This suggests the change of variables 



u - 



4tt 2 

a~ 



The transition thus happens when s 



(4.59) 



which corresponds to the shift of the argument in the theta-functions by imaginary quarter- 
period. Using the standard transformation formulas [18], we rewrite the equations in terms of 
the new variable u: 

g = I - \ — In (0 4 (O)0 2 (iiO + 2 (O)0 4 (iu)) ( 4 - 60 ) 



2 2 du 

6 2 (0)6 4 (iu)6 2 (iu) - 6 2 (0)6 4 (0)6 2 (iu) _ flf (0) - <9|(0) + E 2 
e\{iu) 



(4.61) 



60 2 (O)0 4 (O) 

The left-hand-side of the last equation has an extremum at u = 0, because it is an even function 
of u. This extremum is actually a minimum, and it is also possible to show that this is a global 
minimum on the real axis. The solutions are thus real if the right-hand-side exceeds the value 
of the left-hand side at u = 0, which happens at sufficiently small A. At A = A c , the right-hand 
side crosses the minimum of the left-hand-side, the two solutions for u collide and go off into 
the complex plane. This critical point precisely corresponds to hitting the branch singularity 
of F({p,m) described above. On the other hand, g is equal to 1/2 plus an odd function of u. 
Consequently, 

9M = \, (4-62) 



which finally proves the equivalence of (4.56) and (3.3). The appearance of odd powers of u 
also means that the g(\) becomes complex if continued past A = A c (fig. [4). 



To find the numerical value of the critical coupling, we can either use (4.56) in conjunction 
with (4.32), or set u = in (4.61). The resulting equation can be simplified by using the 



theta-function identities, and takes the form: 



4E 2 (pi)-2E 2 (pl) + E 2 (p c ) = 0. 



(4.63) 
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Numerically, 



Vc~ 0.10765, 



X r » 35.425. 



(4.64) 



Is A c big or small? The expansion parameter of perturbation theory is commonly identified 
with A/(47r 2 ), which at the critical point is « 0.90. This is of order 1. However, the quantities 
that we can calculate with the help of localization appear to have an expansion in powers of 
-8?r 2 /A ^y e ]2 ave identified this expansion with the OPE in the underlying field theory. 



P 

The coefficients of this expansion do not grow too fast (eqs. (4.48), (4.50)). In fact, for the 
free energy the radius of convergence of OPE is p = 1, as follows from (4.47). For the largest 
eigenvalue g(X) the radius of convergence is exactly p c . We have indeed checked that the first 
four-five orders of OPE (eqs. ( |4.48 ), (4.50)) describe the exact functions /(A), g(X) with better 
than a 1% accuracy in the whole allowed range of couplings < A < A c (despite the fact that 
g(X) has a branch singularity at A = A c ). We thus conclude that, while the phase transition 
takes place at A/(47r 2 ) = 0(1), it lies within a perturbative regime in the OPE series. 

A question of interest concerns the order of the phase transition. For that we need to know 
the free energy at A > A c . As we have not solved the matrix model in the strong-coupling 
phase A > A c , we are not able to determine the order of the transition with all the certainty. 
Our numerical results suggest that the transition is rather weak, perhaps of the third order as 
usually happens in matrix models [?]. 

The behavior of g near the critical point also suggests that the transition is rather weak: 
expanding (4.60), (4.61) near u = we find: 



g(X)= 1 --c(X c -X)U. 



(4.65) 



where c is a numerical constant c > 0. 



5 Strong-coupling phase(s) 

Although we do not know how to extend the above analytic solution past the transition, the 
eigenvalue density can be obtained numerically. Figs. [5j |6j [7] show the eigenvalue density at 
different A. We find that fi is a monotonic function of A. The different phases, that emerge 
as A is gradually increased, can be understood in terms of thresholds occurring at \x = nM/2, 
with n = 1,2,3,.... They represent primary and secondary resonances, as explained in section 
3. The picture that emerges is as follows: 

• For A < A c , the eigenvalue density has a shape similar to l/(iryf [i 2 - x 2 ), with the inverse- 
square-root singularities at the boundaries. 

• As A becomes greater than A c , \i overcomes the first threshold, /i c = M/2. The density 
ceases to be a smooth function on the interval and develops two cusp-like features 
at -M + fi and M - fi, which travel towards x = as A is increased. The cusps can be 
understood as resonances due to nearly massless hypermultiplets. Indeed, ±(M-/x) are 
exactly in resonance with the boundaries of the eigenvalue distribution. 
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• The physical reason for the weakness of the phase transition is that at A c the cusps emerge 
from the boundary of the eigenvalue distribution with zero amplitude. Therefore the 
eigenvalue density and observables directly related to it, such as (x 2 ), are not drastically 
modified when crossing from A < A c to A > A c . The amplitude of the cusps - related to 
the number of light eigenvalues - grows as the cusps move towards the interior of the 
eigenvalue distribution with increasing A. 

• The next phase transition happens when \x = M. Numerically we find that the critical 
coupling is A c « 83. At this point, the first two cusps collide at x = and two new cusps 
at x = -fi + 2M and x = fi - 2M are formed near the boundary. 

• As A is further increased, one goes through new phases where more cusps are formed 



(3) 



150, 



pairwise, whenever n crosses nM/2, n = 1, 2, 3, .... The next transition occurs at Ac 
where \i = 3M/2. At the same time, away from the cusps, for large n the eigenvalue density 
approaches the shape of the Wigner semicircle's law, \J p? - x 2 . The amplitudes of the 
cusps also become smaller and smaller as fi grows. 



The Wigner distribution at strong coupling (2.10 ) is the result of averaging over an infinite 
number of infinitely weak cusps. 



Figure 5: Eigenvalue density obtained numerically, a) Subcritical phase with A = 30. b) Supercritical 
phase at M/2 < fj,< M with A = 60. Two cusps formed at ±(M - fi) which move towards x = as A is 
increased. 

In the phase A < A c , the numerics matches the analytic results of the previous section with 
high accuracy. For example, at A = 30, 

(x 2 ) num « 0.11547 , (x 2 ) analytic « 0.11545 , (5.1) 
{x\ um « 0.01878 , (x 4 ) ana i ytic « 0.01876 . (5.2) 

While we could not find p in a closed form, there are quantities which admit simple expressions 
in terms of theta functions. Examples are (x 2n ) discussed before. Another example is 

*»=£(f4). (5-3) 



where £2, £3 were given in (4.32). 
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Figure 6: Eigenvalue density at larger A. Two new cusps form whenever [i crosses nM/2, for any 
positive integer n. a) Phase M < y < 3M/2 with A = 130 and cusps at fj,±M and ±(2M-/i). b) Phase 
3M/2 </j<2M with A = 180 and cusps at n± M, ±(2M - y) and ±(3M - ji). 
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Figure 7: At very large A the eigenvalue density approaches a semicircle distribution superposed with 
multiple cusps whose amplitude vanishes as A -> oo. a) A = 1500. b) A = 5000. 



6 Conclusions 

We have solved for the large-iV master field of the J\f = 2* SYM using supersymmetric local- 
ization [T]. Interestingly, the free energy and the expectation values of large Wilson loops, 
which are calculable by localization, have only non-perturbative, OPE-like terms in their weak- 
coupling expansion. For the free energy, the OPE has a finite radius of convergence. The first 
terms in the OPE can perhaps be reproduced from a more direct field-theory calculations. 

To our surprise, the phase structure of N = 2* theory appears to be very complicated, with 
consecutive phase transitions occurring as A is increased and accumulating towards A = oo. 
Although large A" phase transitions are common in matrix models, such a complicated phase 
structure is quite unusual and to the best of our knowledge has never been encountered before. 
We were able to solve the model analytically only in the most weakly-coupled phase. It would 
be very interesting to find the exact solution in the other phases that would allow, for example, 
to determine the order of the phase transitions. 
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The strong-coupling regime of M = 2* SYM can be successfully confronted [I] with the 
holographic predictions [3] of the Pilch- Warner geometry [2], owing to localization. Although 
the strong-coupling solution appears smooth and simple, the approach to the strong-coupling 
regime is very irregular, with infinitely many phase transitions taking place on the way. The 
asymptotic strong-coupling eigenvalue density, used in [4] for comparison with holography, 
actually has an irregular fuzzy fine structure (fig. [j]). At finite but large A there are 0(y/\) 
resonance cusps whose amplitudes diminish with growing A. The asymptotic density is a result 
of averaging over these resonances. It would be interesting to understand if irregularities of 
the eigenvalue distribution affect the subleading orders in the strong-coupling expansion. Their 
effect may then be visible in quantum string theory on the Pilch- Warner background. 
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A Theta functions 

The four Jacobi theta-functions are defined as 

oo 

0i(*|g) = 2gi £ (-l)'V (n+1) sin(2n + 1),2 

n=0 
oo 

e 2 (z\q) = 2gi £ <f cos (2n + I) z 

n=0 
oo 

e 3 (z\q) = l + 2^<f 2 cos2nz 

n=l 

oo 

6 A (z\q) = 1 + 2^ (-l)V 2 cos2nz. (A.l) 



n=l 



Along with the theta-functions we also use the Eisenstein series: 



= ^Et^ <'^> 

n=l 



r 

°° ^3g2n 



EM 2 ) = l + 240 £f^ (A.3) 
and the Dedeking eta-function: 

oo 

r](q) = q^Y[(l-q n ), (A.4) 

n=i 

which satisfies 

E 2 ( q )=2Aq f ^. (A.5) 
77(g) 
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The Eisenstein series with index larger than two can be expressed in terms of the theta- 
functions, in particular: 

£ 4(g2) = g§(0| g ) + <9§(0[g) + g|(0[g) _ (A 6) 

The standard elliptic integrals are expressed in terms of the theta-functions via the following 
set of equation^] 

q = e~^~ (A.7) 

m " »!(ok) (A ' 8) 

K = |9|(0|«) (A.9) 

£ = 6 «O0 (A ' 10) 

sin ^ = (A . n) 



7 2 



(O|g)04(v|g) 



The inverse of (A. 11) reads 



The transformation properties of the theta-functions, as well as many useful identities are 
listed in |l8j^J We also make use of the Jacobi identity: 

9l(o\ q ) = d$m + di(o,q)- (A.H) 
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